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Segue 1: the best dark matter candidate dwarf galaxy surveyed by MAGIC 
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Abstract: Despite the interest in Dark Matter (DM) searches is currently more focused on underground experiments, a 
signature of DM annihilation/decay in gamma-rays from space would constitute a smoking gun for its identification. In 
this contribution, we start with a brief review of the efforts of the ground-based MAGIC Cherenkov telescopes system 
to detect DM signatures from dwarf satellite galaxies orbiting the Milky Way halo. We then present the recent survey of 
Segue 1, considered by many as possibly the most DM dominated satellite galaxy known in our galaxy. No significant 
gamma-ray emission was found above the background in around 30 hours of observation. This is the largest survey ever 
made on a single dwarf by Cherenkov telescopes. We present a novel analysis that fully takes into account the spectral 
features of the gamma-ray spectrum of specific DM models in a Supersymmetric scenario. We also discuss the prospects 
of detection after the Fermi observation of similar objects at lower energies. 
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1 Introduction 

Over the last century, a great amount of observational ev- 
idence at different scales has pointed to the existence of 
an unknown, invisible mass component that accounts for 
almost 25% of Universe's energy budget fH. Neverthe- 
less, despite strong efforts over the years, Dark Matter is 
still one of the greatest mysteries of modern Physics. No 
experiment so far has been able to detect DM, directly or 
indirectly, and shed light on its nature, although recent di- 
rect detection experiments might have found a marginal 
signature of DM J2j[3]]. Experimental evidence indicates 
that the DM particle should be a weakly interacting mas- 
sive particle (WIMP). WIMPs were thermally produced in 
the early Universe, are stable in time and of non-baryonic 
nature |4|. Through self-annihilation, WIMPs should pro- 
duce Standard Model (SM) particles, and some of the by- 
products, like photons, hadrons and leptons, might be ob- 
servable from Earth. Supersymmetry (SUSY) provides a 
natural candidate for the DM WIMP - the lightest SUSY 
particle neutralino. 

Given the natural mass range for the SUSY WIMPs (be- 
tween a few GeV and a few TeV) Q, it may be possible 
to find, via very energetic photons, some signatures of the 
DM annihilation in the energy range of the Imaging At- 



mospheric Cherenkov Telescopes (IACTs). The typical an- 
nihilation gamma-ray spectrum is predicted to be continu- 
ous and featureless, due to the gammas mainly being pro- 
duced from pion decays and final state radiation of charged 
particles. Nevertheless, some distinctive spectral features 
could be present, like the line emission (from WIMP's an- 
nihilation into a pair of gammas or a gamma and a Z- 
boson), a cut-off or a spectral hardening due to internal 
bremsstrahlung, all dependent on the mass of the DM par- 
ticle. 

The gamma-ray flux from DM annihilation that is expected 
to be measured on Earth is determined by two independent 
factors: one coming from particle physics and the other 
related to the astrophysics. The particle physics term gives 
the number of gammas, produced above the given energy 
threshold for a certain DM model, and it does not depend 
on the source. On the other hand, the astrophysical factor is 
determined by the DM density profile of a specific source, 
as well as its distance. 
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2 Dwarf Spheroidal Galaxies (dSphs) 

Given that the annihilation rate is proportional to the 
squared DM density, best places to look for WIMPs are 
those regions with large concentration of DM. 

Dwarf Spheroidals, the satellite galaxies of Milky Way, 
have been identified as one of the most suitable candidates 
for DM searches: their ratio between the mass inferred 
from gravitational effects and the mass inferred from lu- 
minosity (M /L) can reach values from 100 up to few 1000 
Mq /Lq; they are relatively close and they are expected to 
be free of astrophysical objects whose gamma-ray emission 
might 'hide' the DM signal HQ. 

Until several years ago, only nine dSphs were discovered 
orbiting the Milky Way. These, so called 'classical' dwarfs, 
were joined by 1 1 new, ultra-faint objects thanks to the 
discoveries of the Sloan Digital Sky Survey (SDSS) from 
2005 on J8). Among these object, several have been identi- 
fied as excellent candidates for indirect DM searches with 
MAGIC. 

3 MAGIC Observations 

MAGIC is a ground-based system of two, 17 m diameter 
IACTs, located at the Observatory Roque de los Mucha- 
chos, in the Canary island La Palma (28.8° N, 17.8° W, 
2200 m a.s.l.). MAGIC-I has been operational since 2004, 
and in 2009 it was joined by MAGIC-II. Together, in 
stereoscopic mode, they allow for observations of signif- 
icantly improved sensitivity, lower energy threshold and 
better energy and angular resolution [9 1. 

MAGIC has observed three dSphs so far: Draco [10|, 
Willman 1 [11| and Segue 1 [12|. All of these observa- 
tional campaigns were carried out when only the first tele- 
scope was operating. Neither case resulted in detection 
of gamma-ray flux above the background. The flux up- 
per limits (ULs) were calculated, above a certain energy 
threshold, for each dwarf and used to probe the parameter 
space of mSUGRA. mSUGRA is a 5-dimensional subspace 
of the Minimal Supersymmetric Standard Model (MSSM), 
defined by universal masses of gauginos (jn 1 / 2 ), scalars 
(mo) and trilinear couplings (Aq), as we H as by the Hig- 
gsino mass term sign (p) and by the ratio of the vacuum 
expectation values of two Higgs fields (tan/3) [13|. For 
most mSUGRA models, the lightest supersymmetric parti- 
cle and DM candidate is the neutralino, a linear combina- 
tion of the super-partners of the neutral Higgs bosons and 
gauge bosons. 

Findings of these analyses are discussed, for each dSph 
separately, in the following paragraphs. Segue 1 results are 
discussed in more detail. 

3.1 Draco observations 

Draco is one of the nine 'classical' dSphs, located 82 kpc 
away from us. From the kinematics of over 200 stars as- 



sociated to this galaxy, it has been derived that its M/L is 
around 200 M Q /L . 

Draco was observed with MAGIC-I for 7.8 h in 2007. A 2a 
upper flux limit on the steady emission of l.lxlO -11 ph 
cm -2 s _1 was found, for energies above 140 GeV, assum- 
ing a generic annihilation spectrum without a cut-off and 
with a spectral index of -1.5. For the purpose of constrain- 
ing the mSUGRA parameters space, for several benchmark 
models, we calculated the integral upper limits on expected 
gamma-ray fluxes for energies above 140 GeV and com- 
pared these values to experimental ones. Results showed 
that ULs obtained from MAGIC Draco observations are 
C(10 3 -10 9 ) away from constraining the mSUGRA phase- 
parameter space. 

3.2 Willman 1 observations 

Willman 1 was recently discovered as an ultra-faint dSph, 
located 38 kpc from us, with M/L ~ 500-700 M Q /L G . 
Willman 1 was observed by MAGIC-I in 2008 for 15.5 h. 
No significant gamma-ray emission was found above 100 
GeV, corresponding to 2a integral ULs of order of 10~ n 
ph cm~ 2 s _1 . Assuming the NFW profile of the DM halo, 
flux was calculated for four different, mSUGRA modified 
benchmark models 1 14 1. ULs on the Willman 1 flux mea- 
sured by MAGIC-I were still O(10 3 - 10 5 ) far from putting 
constraints on the considered DM models. 

3.3 Segue 1 observations 

Segue 1 is a ultra-faint object discovered in 2007 by the 
SDSS. Located 23 kpc from the Sun, it was initially con- 
sidered a globular cluster, but recent results from the kine- 
matics of its 66 stars favour a dSph nature ifTBll . In that 
case, Segue 1 is the most DM dominated object known so 
far, with M/L ~ 1320-3400 M /L Q fl6l . 

The Segue 1 observations were carried out with MAGIC-I 
from late November 2008 until March 2009, for a total of 
29.4 hours of good quality data, at low zenith angles, guar- 
anteeing a low energy threshold. This is the longest survey 
on one dSph conducted by IACTs so far. No significant 
gamma-ray signal was found above the background. 

We calculated the integral ULs, for 1 1 energy bins (above 
100 GeV), assuming a power law spectra (with spectral in- 
dex taking values between -1.0 and -2.4). Obtained values 
for the ULs are of order ~ 1 x 10 -11 ph cm _2 s _1 (above 
100 GeV) and - 2 x 10~ 12 ph cm^s- 1 (above 200 GeV). 

3.3.1 Constraints 

In order to compare our Segue 1 results with the predic- 
tions from mSUGRA, we calculated the expected gamma- 
ray flux for different models. 

The astrophysical term of the flux is model independent, 
and to estimate its value, we assumed Einasto profile flTl . 
Taking into account the source's extended emission, our 
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Figure 1: Annihilation cross section ULs from Segue 1 
MAGIC data computed for individual points in the scan. 

integration region enclosed ^64% of the total flux from 
DM annihilation, so the total and effective astrophysical 
factor had values of J(Afi) = 1.78xl0 19 GeV 2 cm" 5 sr 
and J(AO) = 1.14xl0 19 GeV 2 cm~ 5 sr, respectively. 

As for the particle physic term, to fully study the mSUGRA 
phenomenology, we performed a grid scan over the param- 
eter space: for neutralino masses in the detection reach 
of MAGIC-I, and for both positive and negative sign, we 
varied the mSUGRA parameters, linearly and in 40 steps 
over the chosen ranges, to finally obtain over 5 x 10 6 mod- 
els. Each model was then tested with DarkSUSY 5.0.4 
for being unphysical and for its compliance with the SM 
experimental constraints. Also, for every model, using 
Isasugra 7 . 78, we calculated relic density and com- 
pared it to the value derived by WMAP [ 18 1 within 3 times 
its experimental error aw map- 

For each DM model in the scan, we computed integral 
flux UL (above energy threshold Eq), using Segue 1 data 
and the specific gamma-ray spectrum of the particular DM 
model. Since the spectra for the various models have dif- 
ferent shapes and cut-offs, the energy Eo from which the 
integral flux {$>e>e ) was calculated was chosen indepen- 
dently for every tested model. We selected the value of Eo 
which minimized the expected flux upper limit, in the case 
that no excess of events was found. 

In order to put constraints on the DM annihilation rate, we 
translated the calculated ULs (after the Eq optimization) 
into averaged cross section (<r ann v) UL . We then plotted 
those values (Figure 0, as a function of m x , together with 
the model expectations. Each circle on the plot represents 
one model, and for each model, there is one UL (square). 
All the models in the scan for which the relic density is 
smaller than the value from WMAP are plotted as crosses 
for these cases, DM would have to had been produced non- 
thermally, in order to reach correct relic density. 



Figure 2: ENFs as a function of the DM mass models in 
the scan providing a relic density compatible with WMAP 
scan (circles) or below (crosses). 

Due to the dependence of ULs on the energy spectrum, we 
found it useful to provide the exclusion plot in terms of 
enhancement factor (ENF) - a parameter equal to the ra- 
tio between the UL on the averaged cross section and the 
value predicted by mSUGRA. ENF basically indicates how 
much intrinsic flux boost particular model needs in order to 
be detectable. Figure|2]shows the ENF values, as a function 
of m x , for models whose relic density is compatible with 
WMAP value or is below it. ENF distribution for these 
two cases is plotted in the upper panel. It can be seen that 
models with lowest ENF are the ones with highest annihi- 
lation cross section and neutralino mass above 200 GeV, 
while the majority of points has the ENF greater than 10 4 . 
For the models compatible with WMAP the lowest ENF is 
of order 10 3 . It is noticeable that their ENF distribution 
is quite wide this is a consequence of the very high ULs 
relative to the models with small neutralino masses, and 
in general to the models whose gamma-ray flux above the 
energy threshold is low. On the other hand, models with 
relic density below the WMAP value have their annihila- 
tion cross sections closer to their upper limits. 

It is interesting to compare the Segue 1 results of MAGIC- 
I and Fermi/LAT: due to different energy ranges, Fermi is 
better suited for low mass neutralinos, while for the heav- 
ier DM particles searches with MAGIC are preferred. With 
the currently available data, Fermi already excludes some 
mSUGRA models with light neutralino and relic density 
smaller than the one provided by WMAP [19]. With 5 
years of data it is expected that Fermi could probe some 
points studied here, and, in case of no detection, exclude 
most of the large ENF cases. That would leave the space 
of high-mass neutralinos and relatively small ENF values 
unconstrained, and that is precisely the region suitable for 
MAGIC. 
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Figure 3: Exclusion lines for a particular annihilation chan- 
nel and the correspoding regions that provide a good fit to 
PAMELA data. 



For the completeness of this study, it should be pointed out 
that contribution of monochromatic lines to the annihila- 
tion spectrum has been neglected, as expected to be sub- 
dominant compared to the continuous emission. In addi- 
tion, potential presence of substructures in DM halo has 
not been considered. Their effect on the flux increase is 
believed to be by a factor few at most. Finally, the parti- 
cle physics term in the predicted flux does not include the 
Sommerfeld enhancement. This mechanism is expected to 
boost the flux value, in some cases, even by a factor of 10 4 . 

3.3.2 Impact on the PAMELA preferred region 

We have considered the interpretations of the PAMELA 
positron excess |20| in the light of DM. According to some 
DM models, in order to fit PAMELA data, the DM particle 
should be heavy and annihilate into many leptons, mainly 
/i + ^t~ and t + t~. Additional annihilation channel through 
an intermediate state <j>, that decays into e + e~, is also con- 
sidered, although the annihilation cross section required for 
this process is rather high, e>(lCT 23 - 1CT 24 ) |2H . 

We compared our upper limits, obtained for Segue 1 data, 
to the regions in the (m x , (<r ann v)) plane that provide a 
good fit for PAMELA data (Figure |3). Annihilation chan- 
nels /J + /i~, t + t~ and (fxfr (with = 1 GeV) were con- 
sidered. These regions have been adapted from |22|, after 
rescaling from a local DM density (from 0.3 GeV/cm 3 to 
0.43 GeV/cm 3 ) (23 1. It can be seen that ENF needed for 
our ULs to probe the PAMELA-favoured regions is much 
smaller than in the case of mSUGRA. Furthermore, in the 
case of t + t~ annihilation channel, our ULs are probing 
the relevant region. However, considering the 2<r uncer- 
tainty in the astrophysical factor, we do not claim exclu- 
sion of the of the PAMELA region with our Segue 1 data. 
Perhaps in the future, with more more precise value of J, 



we will be able to place that constraint, at least for massive 
DM candidates annihilating in t + t~. 

4 Conclusions 

The observations of dSphs carried out by MAGIC-I have 
not resulted in detections. The acquired data were used 
to calculate the integral upper limits on the annihilation 
gamma-ray flux, and then to test (different) mSUGRA 
models. In the most optimistic case, for Segue 1 data, 
the minimal boost of the DM flux needed for detection is 
ENF > 600, for energies larger than 100 GeV and for mod- 
els whose relic density is compatible with one derived by 
WMAP 

Segue 1 ULs have also been used to probe DM models 
which try to explain PAMELA data. Current results are 
probing the PAMELA preferred region for the case of DM 
annihilation into t + t~, but solid exclusion might be con- 
firmed only if the uncertainties in the astrophysical factor 
for Segue 1 are eventually reduced. 

Although the current observations did not find hints of 
a signal, and are far from putting the constrains on the 
mSUGRA parameter space, we expect improvement for 
observations of dSphs carried out with MAGIC in stereo- 
scopic mode. 
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